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superconductivity We show that anisotropic s-wave gap has the same 

low voltage power law conductance and two peak structure in the density 

of states as d-wave superconductors. The assymetric tunneling conduc- 
es 

tance and zero bias conductance for the c-axis tunneling is shown to occur 
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5_i ' because of finite band splitting coming from the interlayer hopping param- 
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Anomolous single particle tunneling characteristics in NIS and SIS junctions of the 
high T c superconductors have remained a subject of great interest. The single particle 
tunneling conductance in both normal and superconducting states is a measure of the 
density of states in the normal and superconducting states, and so in principle one gets 
a lot of informations about the details of the superconducting gap parameter. This is 
a matter of great current interest, in view of the recent controversy about d-wave or 
anisotropic s-wave symmetry of the gap function. 

Whereas some experimental results ( absence of Hebel-Slichter anomaly at T c |1|, 
Raman intensity and low temperature penetration depth measurements 0) are in 
favour of a d-wave symmetry gap function, the recent experiments (absence of Josephson 
current [§]], temperature dependent gap anisotropy ||, a.c conductivity measurements 
H) are in favour of an anisotropic s-wave gap function. 

It has been argued that JF|, the tunneling characteristics, specifically the quadratic 
rise of current with voltage at low enough voltages and the two peak structure seen in the 
conductance voltage characteristics is an evidence for a d-wave superconductor. In view 
of this, we investigated the single particle tunneling characteristics for the anisotropic 
s-wave superconductors , recently proposed by Chakraborty et. al. || 

Before we come to the specific problem we address, we highlight the main puzzling 
features seen in the tunneling spectroscopy of the high T c superconductors |J . 

(1) At low temperature, V — tunneling current is zero for tunneling along the ah 
axis, but nonzero along c axis tunneling. (2) The ab plane conductance becomes smooth 



at larger temperatures but the c axis conductance goes on incresing with temperature. 

(3) At low temperature and large bias, ab plane tunneling conductance decreases or 
saturates whereas the c axis conductance goes on increasing roughly linearly with voltage. 

(4) The ab plane tunneling shows conventional gap like structures but the c axis tunneling 
does not show any such charecteristics. The c axis I-V characteristics shows a much 
broadened shoulder at the gap edge. (5) Both NIS and SIS junction shows assymetric 
I-V characteristics with respect to the sign of the bias voltage. (6) Both direction 
tunneling shows finite density of states for V < A even at the lowest temperatures. 
For very low voltage current has a quadratic rise. (7) There is large broadening of 
gap in voltage , and a conductance overshoot for c axis tunneling. Trying to explain 
this gap broadening due to inelastic scattering leading to quasiparticle lifetime effect 
can explain the gap broadening, but then the zero bias conductance comes out to be 
much larger than observed. Invoking strong coupling corrections does not help, for 
even though gap broadening of correct magnitude can be obtained, but it is difficult to 
explain the conductance overshoot near the gap edge. (8) For ab axis tunneling, zero bias 
conductance is zero. So there is no density of states at the Fermi energy, but for very 
small V < A there is finite current, showing that there is no fully developed gap, or the 
gap is highly anisotropic with A^ being very small in a substantial region of the Brillouin 
zone. For the c axis tunneling the most common explanation for the ubiquitous zero bias 
conductance and the characteristic V shaped conductance versus voltage characteristics 
is explained , as due to either because of tunneling through localized states in the barrier 



or due to scattering by magnetic impurity inside the junctions. The first process is known 
to give rise to non trivial energy dependence of the tunneling probability and can lead 
to zero bias conductance, and the second process was invoked by many people to explain 
non linear current voltage characteristics in tunnel junctions above the gap. It is worth 
emphasizing that these peculiar features of the c axis tunneling are seen in point contact 
and break junction measurements also. 

In this paper we emphasize on the distinction between the ab plane and c axis single 
particle tunneling channels for both supeconductor to normal and superconductor to 
superconductor (NIS and SIS) junctions. Specifically we shall consider a layer material 
like YBCO or Bi-2212 material. Generalisation to multilayer system is trivial. We model 
such superconductors by two planar BCS superconductors coupled by a single particle 
hopping term along the c axis. We consider also the case, when over and above the single 
particle hopping, there is a Josephson coupling between the planes. We propose that the 
observed assymetry of the normal state in plane tunneling conductance, with respect to 
the sign of the bias voltage in MIC (metal- insulator-cuprate) junction is a consequence 
of the existence of nonbonding and bonding band with finite splitting between them, in 
the cuprates. So far there is no agreement on this observed assymetry. Barrier shape 
effects cannot explain it, because it requires an unusually low barrier height. In the split 
band picture, when the metal is held at positive bias with respect to the cuprate, then 
there are two channels of elastic tunneling into the nonbonding and the bonding bands. 
In the reverse bias situation, only one of the bands takes part in tunneling. So the 



conductance will be assymetric, for MIC and NIS(here both sides are cuprates, but one 
is superconducting and the other is in the normal state) junctions. On the other hand 
for SIS and CIC junctions the conductance voltage characteristics will be symmetric. 

Now for a CIC , NIS or for a SIS junction, when both sides of the junctions are 
cuprates, there is an important difference between tunneling along c and ab axis. In the 
ab axis tunneling geometry electrons tunnel only from antibonding to antibonding and 
bonding to bonding bands. Whereas for the c axis tunneling there is another additional 
channel for conduction, i.e from antibonding to bonding band. This tunneling will be 
present even in absence of a finite bias voltage either way. The chemical potentials for 
the two bands differ by 2t±, where t± is the c axis hopping amplitude. So the tunneling 
along c axis will show a zero bias conductance, but the ab axis tunneling will have zero 
conductance at V = and T = 0. 

We find that the zero bias tunneling conductance observed along the c axis tunneling 
increases with temperature and do not show any sign of saturation at all. This is our 
main result. 

We shall also discuss, the reason why for T > T c the ab axis tunneling characteristics 
becomes smooth, while the c axis tunneling continues to be temperature dependent and 
rises with temperature. Lastly we predict that for MIC geometry tunneling (below T c ) 
the assymetry (or alternatively the background conductance) will be more for lesser 
value of the gap in the superconductor. 

To start we take the effective hamiltonian proposed by Chakraborty et. al. M 
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In this model, there is no hopping term along the c axis from plane to plane, even 

though the band theory estimates for the c axis hopping amplitude t± is about | to 

1 of the inplane hopping parameter. The reason is supposed to be that, due to strong 

correlation in the plane itself, the single particle band motion is absent in the c direction. 

The conduction along c axis is purely due to incoherent processes. This phenomenon is 

termed as "confinement" by Anderson. The net effect of being, that there is no need 

to keep the single particle hopping term along c axis. On the other hand coherent 

propagation of "singlet objects" ( pairs of elctrons)is possible. That is the origin of the 

last term ( Josephson coupling of a very unusual kind). It should be emphasized that, 

it has not been proved within a realistic model for hig T c superconductors. 

We prefer to keep the band term in the hamiltonian. The origin of subbands can 

be understood as follows. We consider a two layer material like YBCO. Th individual 

layers can be modelled by a 2-d tight binding band with dispersion, 

€k = — 2t(cos(kx) + cos(ky)) + 4t'cos(kx)cos(ky) 

where t and t' are nearest and next to nearest neighbour hopping in the planes of some 
effective site. We take, t = 0.3eV and 4t± = 0.45eV. For two closely spaced planes, in 
interlayer matrix element 

t±{k) = t±(cos(kx) — cos(ky)) 2 



results in formation of subbands, 



E M {k) = e(k)±t ± (k) 



Where the (f> and ip are antibonding and bonding band fermions defined as, 



The location of chemical potentials will be detrmined by the doping. 

To illustrate the difference between tunneling along ab axis and along c axis we 
write down the tunneling hamiltonian without any explicit dependence of the tunneling 
amplitude on momenta or energy For tunneling along ab axis, the hamiltonoan will be 

kp 

where a denotes layer index (1 and 2). The c axis tunneling hamiltonian on the other 
hand will be 

J2 k P T k P clA, + 1^2= j2 T k P (4,<t>ka +<p->i> + £^(4^ + h.c) 

kp kp 

It is clear, that for c axis tunneling there is an extra channel for conduction , i.e from 
the nonbonding to bonding band which is absent for the ab axis tunneling. 

MIC junction tunneling . 

To illustrate the effect of this band splitting in the normal state itself , let us consider 
the MIC junction tunneling. Within the independent electron approximation, the single 
particle tunneling current is given by, 

/ = [\T\ 2 N 1 (E)N 2 (E + eV)[f(E)-f(E + eV)]dE 



where energy E is measured from chemical potential, iVi and N2 are density of states of 
two electrodes (one usual metal and the other one being the cuprate ), V is voltage bias 
and / is the Fermi function. If we take the density of states of the metal ri2=const and 
that of the cuprate to be Ni(E) = N^ + N 1 p8(—2t± — E) where N^ and N^ can be taken 
as constants for simplicitythen from the above two equations , we get, 

with a = N^/Nq. For T = 0, G(V) = 1 + a. It is clear that, the tunneling conductance 
is assymetric with respect to bias. For positive bias, G(V) ~ {T^N^l + A) and for 
negative bias, « |T| 2 iV^ if \V\ > 2t± and « |V| 2 A^(1 + A) otherwise. This is true at 
T = 0. At nonzero temperatures , the conductivity assymetry will be seen at lower bias 
and the absolute value of conductivity will decrease. 

At this point, we compare our model with that of Levin and Quader |ID |, who 



also consider a split band picture. We insist that there is a major difference between 
our viewpoints as regards the role of the split bands. Levin et. al. 1(J assume that the 



bonding band(^ band) is almost submerged below the Fermi surface. For the underdoped 
case, the ip hole band is completely filled and frozen much below the Fermi surface, and 
do not take part in tunneling to the metal on the other side of the junction. Consequently 
there will not be any conductivity assymetry for underdoped case. For larger doping case, 
both the bands will be partially filled and take part in tunneling. Moreover one needs 
additional assumptions, that the ip band is actually a band of nondegenerate band of 
fermions since their number is so small. One needs to have, in an adhoc fashion, different 



dispersion for and ip fermions(linear and quadratic in momenta) to reproduce some 
normal state properties. This picture is approximately right when t± is large, giving rise 
to large band splitting. We assume , on the other hand that the band splitting is small 
(small t±). So, even at small doping concetrations , both bands will be partially filled. 

Within the interlayer tunneling mechanism of superconductivity, even though the 
intralayer BCS coupling gives a small T c = 5K on its own, a very small t± is enough to 
raise the T c to large values = 90K through the Josephson coupling term. We have not 
yet made a detailed study of the doping dependence on the MIC tunneling. In other 
words a small band splitting explains the observed assymetry in tunneling conductance 
at small doping as well as very high T c in these materials. 

CIC junction tunneling 

For CIC junctions when both the electrodes are high T c materails (break junctions), 
a look at the tunneling hamiltonians for the ab and c axis tunneling shows that, for ab 
tunneling, the elctrons tunnel from to and from if) to if) bands only. For the c axis 
tunneling, cross tunneling also takes place. If T^ T^ and T^ are the tuneling matrix 
elements between the respective subbands of both electrodes, we get 

G ab (V) = \T H \ 2 NH(l + aj) - *lf C 2t ^~ V )] 
and 



G C (V) = \T„\*N}[(1 + a\ + 2a 2 ) ■ - a 2 /( ^ ) - (a 2 + a{) /(- ^ n 

where a\ = N^T^/N^T^ and a 2 = N^T^/N^T^. The main features of this expres- 
sion are: (1) The conductance voltage characteristics is symmetric with respect to bias 



for both ab and c axis tunneling. (2) For c axis tunneling, there is a zero bias current 
coming from cross tunneling, which is operative even at zero bias because of finite band 
splitting. For ab axis tunneling there is no zero bias current. (3) There is a zero bias con- 
ductance for both ab and c axis tunneling At T = and V — 0, G(V) = T%,Ni(l + af). 

If we take the tunneling matrix element \T\ 2 = 1 + (V/V c ) 2 , then the conductivity 
increses with voltage, but with different slopes for positive and negative bias for the 
MIC junctions and with same slope for CIC junctions. The 1 + (V/V c ) 2 dependence 
of T 2 comes because of Coulomb blockade effects in the junctions. One can also get 
a linear conductance for small voltages due to inelastic scattering in the junctions as 
we mentioned earlier. These kind of approaches are specially tailor made to explain 
the linear conductance in the cuprates. As we pointed out that the linear conductance 
is observed only for the caxis tunneling, one has to explain why, inelastic scattering 
and coulomb blockade effects are not seen for the ab axis tunneling also. Moreover the 
ubiquitous linear conductance is seen in point contact tunneling also. We do not attempt 
to explain this important feature here. The main thrust of our arguments is to show the 
natural origin of the tunneling assymetry in the high T c materials. One more important 
consequence of our model is that in the superconducting state, there is a finite zero bias 
conductance for c axis tunneling. For the inplane tunneling this is absent. This will be 
explored next in SIS and NIS junction tunneling geometries. 

SIS and NIS junction tunneling 

The mean field hamiltonian in the superconducting phase is, 



(2) 



The hamiltonian looks like a sum of two BCS reduced hamiltonins for the bonding and 
antibonding electron systems. The generalised gap equation will be 

1 l^tanh((3E$/2) 1 _ tanh(l3Ef/2) 



(V+'-f) *r 2Ef 2^ 2Et 

where, 

E£* = y/(e k ±t ± )* + A* 

Going from summation to integral and converting to energy variables it is not very 
difficult to see that the T c is given by 

k B T c = JujZ-tl—e'ww+'lM (4) 

IT 

for small values of t±, where e 7 = 1.781. 

We have solved the gap equation numerically for different temperatures. In NIS 
junctions, the tunneling current is given by, 



I NIS = = -£\ T \ 2 [ulS(eV+E k -^)[f(E k )-f^ p )] + v l5(eV-E k -Q[l-f(E k )-f(Q}} 

kp 

For the SIS junction the corresponding expression is, 



Isis = Y.\ T \\ l -f( E k)-f{E P )){vlul8{eV-E p -E k )-ulvp{eV + E p + E k )) 

kp 

+ (f( k ) - f(E p ))(u 2 k u 2 5(eV + E k - E p ) - vk 2 v 2 p 8(eV + E p - E k ))\ (5) 
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The normalised conuctance versus voltage for the ab and c axis tunneling are plotted in 
Fig.l and Fig. 2 respectively. The notable features are, (1) At T = there is a sharp 
voltage threshold for conductivity for the ab axis tunneling, whereas there is a finite 
zero bias conductance for the c axis tunneling. (2) The sharp voltage threshold for ab 
axis tunneling gets washed out at a very small temperature (4-ff). (3) The plots for 
conductance at T = 4, 10 clearly shows the characteristic two peak structures seen in 
experiments. For d-wave superconductors also one gets similar two peak structures. 

In Fig. 3 we plotted the current versus temperature for ab axis tunneling for 2 and 
5 degree Kelvin. We emphasize that, even at very low temperatures ( 5 K) the current 
rises quadratically with voltage at very low voltages. This clearly shows that, the gap 
at most places of the Brilloiun zone is very small and falls faster with temperatures, 
than usual BCS temperature dependence of gap. Thus at small but finite temperatures 
the anisotropic s-wave superconductor becomes indistinguishable from a superconductor 
with gap nodes. 

Fig. 4 shows the temperature dependence of the normalised zero bias conductance 
for both ab and c axis tunneling. One extraordinary feature of the interlayer tunneling 
gap function is that the gap along V — M direction is large and almost temperture 
independent and retains its full gap value at T = up to about 90% of the T c , and 
then falls almost like a weak first order transition. On the other hand gap in any 
other direction falls much faster that the usual BCS gap suppresion due to thermal 
fluctuations. The momenta averaged gap also falls very slowly with temperatures, as 
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observed in the recent photoemission experiments. This is true when Tj > Vj, cs , or in 
other words, when interlayer tunneling is dominant for very strongly coupled layers. 

For weaker Tj or with larger in plane Vb cs , the averaged gap falls faster with tempera- 
tures and slowly approach the usual BCS temperature dependence. Notice that all these 
pecularities are only because of the 1 — k summation in the interlayer josephson coupling 
term, as emphasized by Anderson. For a more conventional josephson coupling, where 
the individual momenta of the partners of the cooper pairs are not conserved, and only 
the center of mass momenta is conserved, i.e a josephson term with double momenta 
summation, we do not get the above mentioned features at all. 

Two things follows automatically from above discussion. One is that, in the interlayer 
mechanism, the gap magnitude in most part of the BZ is very low (1-3 meV) and also 
very fragile as far as thermal fluctuation is concerned. The gap in these regions falls 
faster than in the usual BCS gap. This would mean that we shall not get any sharp 
gap features at all at any finite temperatures in tunneling experiments. This is what is 
observed in our numerical calculations at finite temperatures. In Fig. 3. we show the 
I — V characteristics at T = and T = 5 degrees for tunneling along the ab plane. We 
see clearly that already at T = 5 degrees there is finite current at very low voltages. In 
other words, indeed it will be very difficult to distinguish between, the situation where 
the gap function has gap nodes on the Fermi surface like in d-wave superconductors, 
and the interlayer case. 

Note that for tunneling along, the c axis there will be a finite current for arbitrar- 
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ily small voltages. So in ceramic materials, where we measure some average current 
along both directions, we shall always get an I-V characteristics looking just like a su- 
perconductors with gap nodes on the Fermi surface even at T = 0. For single crystal 
measurements , and for ah plane tunneling there will be a sharp voltage threshold, but 
no sharp threshold for small but finite temperatures. It is worth emphasizing that, we 
do not really know how impurities and inhomogeneities suppress the interlayer tunneling 
gap. 
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Figure Captions 



1. Conductance ( Normalised with respect to normal state) vs Voltage for ab axis 
tunneling for Temperatures 2 (Solid line), 4 (dashed line) and 10 (dotted line) 
degrees. 

2. Conductance ( Normalised with respect to normal state) vs Voltage for c axis 
tunneling for Temperatures 2 (Solid line), 4 (dashed line) and 10 (dotted line) 
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degrees. 

3. Current Vs Voltage for the ab axis tunneling, for temperatures 2 (solid line) and 
5 (dashed line) degrees. 

4. Zero Bias Conductance vs Temperature for ab axis (solid line) and for c axis 
(dashed line). 
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